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Abstract 

This paper deals with bounds satisfied by the effective non-symmetric conductivity 
of three-dimensional composites in the presence of a strong magnetic field. On the one 
hand, it is shown that for general composites the antisymmetric part of the effective 
conductivity cannot be bounded solely in terms of the antisymmetric part of the local 
conductivity, contrary to the columnar case studied in [15]. So, a suitable rank-two 
laminate the conductivity of which has a bounded antisymmetric part together with a 
high-contrast symmetric part, may generate an arbitrarily large antisymmetric part of 
the effective conductivity. On the other hand, bounds are provided which show that the 
antisymmetric part of the effective conductivity must go to zero if the upper bound on 
the antisymmetric part of the local conductivity goes to zero, and the symmetric part 
of the local conductivity remains bounded below and above. Elementary bounds on the 
effective moduli are derived assuming the local conductivity and effective conductivity 
have transverse isotropy in the plane orthogonal to the magnetic field. New Hashin- 
Shtrikman type bounds for two-phase three-dimensional composites with a non-symmetric 
conductivity are provided under geometric isotropy of the microstructure. The derivation 
of the bounds is based on a particular variational principle symmetrizing the problem, 
and the use of y-tensors involving the averages of the fields in each phase. 

Keywords: bounds, homogenization, magneto-transport, strong field, multiphase composites. 

AMS classification: 35B27, 74Q20 



1 Introduction 

It is known since the seminal discovery of Hall in the end of the 19th century [19] , that a low 
magnetic field perturbs the matrix resistivity (or equivalently the conductivity) of a conductor 
by inducing a small non-symmetric part characterized by the so-called Hall coefficient. In the 
80's Bergman [5] first gave a general formula for the effective Hall coefficient involving currents 
that solve the symmetric conductivity equations in the absence of a magnetic field. However, 
there are few explicit formulas for the effective Hall coefficient except in very particular cases 
like two-phase two-dimensional composites [24, 6, 16], or columnar composites [7, 9, 31, 22, 23]. 
The situation is still less favorable in the strong field case [8, 10], namely when the symmetric 
part and the antisymmetric part of the conductivity are of the same order. In three dimensions. 
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only when the antisymmetric part is constant do we have an exact formula for the antisymmetric 
part of the effective tensor [32] . So, rather than trying to get exphcit relations for the effective 
tensors it seems more practical to derive bounds. The theory of the bounds in homogenization 
has gone through a considerable development since the original variational approach of Hashin 
and Shtrikman [20]. We refer to [26] for a comprehensive survey. In fact, very little is known 
about bounds for strong field magneto-transport. Recently, we derived in [15] optimal bounds 
for multiphase columnar composites. The aim of this paper is to extend, at least partially, the 
result of [15] to two-phase three-dimensional composites. 

In the present context we consider a three-dimensional conductor having a periodic structure 
(this is actually not a restrictive assumption) in the presence of a fixed vertical strong magnetic 
field. Under the transverse isotropy assumption the local conductivity of the conductor takes 
the general form 

(a{y) -c{y) \ 
<y) = cijj) aijj) , for t/ e M^ (1.1) 
V h{y)) 

where the coefficient c{y) is induced by the presence of the magnetic field parallel to the 1/3-axis, 
which also infiuences aiy) and h{y) and causes them to be non-equal in the case of a conductor 
that is isotropic in the absence of the magnetic field. Similarly, assuming an transversely 
isotropic microstructure, or at least one that is invariant under 90° or 120° rotations about the 
y3-axis, the constant effective conductivity of the composite is given by 

(a* — c* 0\ 
c* a* . (1.2) 
k) 

Our goal is to derive bounds for the effective coefficients a*, 6*, c* of cr* in terms of the coefficients 
a(y), h{y), c{y) of the local conductivity (T{y). 

In Section 2, we derive elementary bounds (see Theorem 2.1) on the effective coefficients a*, 
6* and c*. These are obtained by taking uniform trial fields in a variational principle for non- 
symmetric tensors deduced in [25, 18] from a symmetrization of the constitutive law j = ere, 
and its adjoint, adapted from the variational approach performed in [17] for complex tensors. 

In Section 3, we show (see Theorem 3.1) that contrary to the columnar case of [15], it is 
not possible to bound the antisymmetric part of the effective conductivity cr* only in terms of 
the coefficients c{y). Indeed, when a{y) is independent of y^ for a vertical columnar structure, 
the key ingredient for the derivation of the optimal bounds in [15] is based on the positivity of 
the (2 X 2) determinant of the local electric field [1, 2], i.e. 

Ai^2{DU) := ^1^1 d2U2 - diU2 d2Ui > a.e. in M^ (1.3) 

where the vector- valued potential U — {ui,U2,us) solves the conductivity problem 

f div (aDU) = in 

\ U{y) — y is F-periodic. 

Due to a suitably constructed rank-two laminate with high-contrast conductivity, we prove 
simultaneously that the inequality (1.3) does not hold, and that arbitrarily large effective 
coefficients c* can be obtained while the local coefficient c(y) is bounded. This negative result 
agrees with the pathologies obtained in [13, 14] with different microstructures, related to bounds 
on the effective Hall coefficient in the low magnetic field regime. As a consequence, a bound 
for c* involves both the upper bound for |c(y)| and the bounds from below and above for a{y) 
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in (1.1) (see Theorem 3.1 and Remark 3.2). This bound shows c* — > when the upper bound 
on \c{y)\ goes to zero, provided a{y) remains bounded from below and above. 

In Section 4, to improve the previous bounds we restrict ourselves to a two-phase local 
conductivity 

(ai -ci 0\ /a2 -C2 0\ 

ci ai + X2{y) C2 a2 , for y e M^ (1.5) 
hj \0 62/ 

with prescribed volume fraction /j = {xi), for i = 1,2, with /i + /2 = 1. We then derive 
(see Theorem 4.1) Hashin-Shtrikman type bounds for the effective conductivity cr*, involving 
three intermediate coefficients ay, by, cy which are explicitly expressed in terms of the entries 
of (Ti, for i = 1, 2, *. In particular, it is shown that the point (ay, —cy) belongs to a disk 
which is tangent to the axis a = at some point (0,c), and which contains the disk passing 
through the points (ai, Ci), (02, C2), and tangent to the axis a = at the same point (0, c) (see 
Figure 1 below) . The derivation of these new bounds is based on a combination of three main 
ingredients: 

• the geometric isotropy of the phases defined in [34] for random composites, 

• the variational principle for non-symmetric tensors, 

• the use of F-tensors similar to [21] (see also [11]), giving relations between the averages 
of the fields in each phase. 

Notations 

• (ci, 62, 63) denotes the canonic basis of M^. 

• / denotes the unit matrix of R^^^, and J := (1 0^ ). 

\ / 

• For any matrix M e M^''^ Af^ denotes the transpose of M, := | (M + M^) the 
symmetric part of M, and M"^ := ^ (M — M^) the antisymmetric part of M. 

• Y denotes the unit cube [0, 1]"^, and (■) the F-average. 

• For a function / defined on the unit sphere S"^ of M^, (7)52 denotes the average of / over 

i.e. 

f 1 r^TT PIT 

{f)s^'—T fiOd^=~r~ / ^'P f {sin 9 cos (p, sin 9 sin (f, cos 9) sin 9 d9. (1.6) 

• For a, 13 > 0, Mj(a,/3;y) denotes the set of the ^-periodic invertible matrix- valued 
functions A : -> R^^^ such that 

VeeR3, A(t/)e-e>«ier and ^"^(y) C • C > ICI' a.e. y e F. (1.7) 

• L'^{Y) denotes the space of the F-periodic functions, which are square integrable in Y. 

• H^{Y) denotes the space of the F-periodic functions, with gradient in Lj{Y)^. 

• Foru-.R'^R, V«:= (fe)i<,<3- 

. For [/ : R3 ^ R3, U = {u,,u,,Us), DU (£t)i<,,<3- 



• For E : R3 ^ R^xs, div (E) := (^) 



l<j<3- 
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2 Elementary bounds on magneto-transport 

To derive bounds on composites we may assume that the associated microstructures are Y- 
periodic (see, e.g., [3] Theorem 1.3.23), where Y is any cube of M^, say Y = [0,1]"^. In this 
section and the next we consider a three-dimensional F-periodic conductor in the presence of 
a strong magnetic field parallel to the y^-Bjds so that the resulting matrix- valued conductivity 
a{y) is given by 

/a{y) -c{y) \ 
^(y) = c{y) a{y) , for y e M^ (2.1) 
V b{y)J 

where the coefficients a{y),b{y),c{y) satisfy for prescribed positive numbers a, a, c > 0, with 
a < a, the following bounds 

a<a{y)<a and \c{y)\<c, a.e. y e M^. (2.2) 

By virtue of the periodic homogenization formula (see, e.g., [4]) the effective conductivity cr* 
associated with cr(y) is given by 

a^: = {aDU), where the potential ?7 solves <| 'lJ[y)^^^ ^ is y periodic i'^-^) 

Recall that cr* is also the homogenized conductivity obtained from the oscillating sequence (7(|) 
as e — )■ by a homogenization process (see, e.g., [4]). 

Now consider a periodic electric field e G L^(y)^ and a periodic current field j £ L^(y)^ 
that solves the conductivity equations 

j — (76, div j — 0, curl e = (2.4) 

and another periodic electric field e' e L^(F)^ and another periodic current field / e L^(y)^ 
that solves the adjoint equations 

j' = a'^e', div f = 0, curl e' = 0. (2.5) 

The average fields are related by the effective tensor cr*: 

0-) = ^*(e), {j')-<7l{e'). (2.6) 



Define the symmetric tensor 



Then, an easy computation yields 



. , eA '■= - [e — e 

es \ _ r 3s \ 



JAJ \eA 




L - =L£;. where ^ ^ (2.8) 



is := 2 (i +/) 



Moreover, mimicking the approach of [17] for complex tensors, extended in [25, 18] (see also 
[26], p. 277) for real but non-symmetric tensors, the following variational principle holds 



ca e Ll{Yf, curl(eA) = 0, (e^) = eo 
js^mYf, div (is) = 0, 05)=Jo- 



, (2.9) 



with the symmetric effective tensor 



[a. 



S\-l 



(^^{(^^)~^ erf — (o"f o-:^^ 



(2.10) 



By substituting constant trial fields = Cq and js = jo in the variational principle one 
immediately obtains the elementary bound 

L, < (L). (2.11) 

This elementary bound implies the following theorem: 

Theorem 2.1. Assuming o"* and (j{y) have the forms (1-2) and (2.1) the constant must 
satisfy the arithmetic and harmonic mean bounds 

1/(1/6) < 6* < {b), (2.12) 

and the pair (a*,c*) must satisfy the circle bounds (which confine (a*,c*) to lie within a circle 
in the a*-c* plane) given by 

(c* - cl^ < (a* - ai) {di - a*) , (2-13) 
where ^ 

aL:={-) , CL:^(-)aL, dL := ( a + -) - ^. (2.14) 



Remark 2.2. Taking the minimum in (2.9) over all fields and js with (e^) = eo and 
(is) = joi ignoring the differential constraints that curl {ca) = and div [js) = gives 
the elementary bound < {L~^). However this does not yield any new inequalities beyond 
(2.12) and (2.13) due to the structure of the matrices and L{y). 

Proof of Theorem 2.1. The proof follows the proof of the elementary bounds in Proposi- 
tion 3.1 of [15]. Assuming cr^, and a{y) have the forms (1.2) and (2.1) we obtain 
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(2.15) 



(2.16) 



The matrix (L) — L* will then be positive semi-definite if and only if (2.12) holds and 

a, > ai, (2.17) 



' <(---)(d, + ^-a.-^) (2.18) 

By multiplying the last inequality by a^a^ and expanding (and using the fact that a^a^ > 0) we 
get (2.13). Also the inequality (2.17) is superfluous as it is implied by (2.13) and the inequality 
dL > aL- 



3 Bounds on magneto- transport: 2d versus 3d 

From (2.2) and the non-negativity of the 4th (or 5th) diagonal element of (L) — we deduce 
the additional (superfluous) bound 

2 / 2 \ —2 

2|c,| <a, + ^ < (a + -) <a + -, (3.1) 
\ a / a 

where to obtain the flrst inequality we have used the fact that x +l/x > 2 for all a; > 0. Thus, 
we have obtained an upper bound on |c*|, but one that involves not only c but also a and a. 
This is in contrast to the case for a columnar conductivity a{y) which is independent of the y^- 
variable, where using the positivity of the determinant of the electric field DU{y) = DU{yi,y2) 
established by Alessandrini and Nesi [1, 2], we proved in [15] that c* satisfies the same bound 
c as the local coefficient c(y) (i.e., |c*| < c). 

Let us now relax the assumption that the effective tensor takes the form (1.2). The composite 
is said to be partially isotropic if the antisymmetric part of cj* satisfies 



(a,)^ = c* J, where J := 1 . (3.2) 




Given a partially isotropic composite we can always subdivide it into square columns with edges 
parallel to the ya-axis and with side length much larger than the existing microstructure, and 
then rotate each square column about its center axis by either 0°, 90°, 180° or 270° with equal 
probability in an uncorrellated way. The resulting polycrystal is invariant under rotations of 
90° about the 1/3-axis and thus will have an effective tensor of the form (1.2) and by a lemma 
of Stroud and Bergman [32] will have the same constant c* as the original partially isotropic 
composite. 

The question naturally arises as to whether for partially isotropic composites |c*| can be 
bounded solely in terms of c, like in the case of a columnar conductivity cr{y)7 The answer is 
no, it cannot. Indeed, we have the following result: 

Theorem 3.1. Consider a periodic conductivity cr(y) given by (2.1) which satisfies the hounds 
(2.2). Assume that the composite is partially isotropic in the sense of (3.2). Then, the effective 
coefficient c* satisfies 

|c*| < - (a^ei ■ ei)^^^ ((7*62 ■ 62)^^^ . (3.3) 
a 

On the other hand, given any arbitrarily large constant /t > there exist transversely isotropic 
conductivities (Jq^^{ii) depending on a parameter 9 > 0, with c{y) G {0, 1} a.e. y E M.^, such 
that as 9 ^ the effective conductivity al ^ is partially isotropic and satisfies 

\im(an)'^ = -KJ or Urn (an = k J. (3.4) 
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Remark 3.2. In the case when cr* is transversely isotropic, taking the form (1.2), the bound 
(3.3) reduces to 

|c*| < — c, (3.5) 
Q. 

and using (3.1) we obtain that 



c / / cl\ C ( (? 



|c.|<- (a + -)--^ <- a + - . (3.6) 
a\\ a I a^j a\ a J 

In contrast to the bounds (3.1) and (2.13) this new bound shows that c* necessarily goes to 
zero as c goes to zero if a and a are held fixed. Also if we add the antisymmetric matrix cqJ to 
(7(y) then the effective tensor will change to cr* + Cq J, implying from (3.5) that the inequality 

|c* + Col < — max (|c+ + co|, |c_ + co|) (3.7) 



a 



holds for all constants cq, where 



c+:=supc(y), c_ := inf c(y). (3.8) 

y&Y 

Taking the optimum value cq = — (c+ + c_)/2 gives the bounds 

12 c* — c+ — c_| < — (c+ — c_). (3.9) 
a 

Remark 3.3. Theorem 3.1 proves that contrary to the columnar case of [15] we cannot expect 
to bound the effective coefficient only in terms of the bound c of the local coefficient c{y). 
Actually, (3.4) shows that arbitrarily large (positive or negative) effective coefficients c* can be 
derived although the local coefficient c{y) only takes values in {0, 1}. Here, the contrast of the 
symmetric part of the conductivity plays a crucial role as suggested in the bound (3.3). This 
is strongly linked to the fact that the (2x2) determinant 

Ai,2(-Dt/) := diUid2U2 - (?iM2(?2Mi, for U = (mi,M2,«3), (3.10) 

does not always have a constant sign throughout the material (see the proof of Theorem 3.1 
below) contrary to the columnar case. 

Proof of Theorem 3.1. 

Bound for c*; The div-curl lemma of Murat-Tartar (see [33, 29, 30]) and the formula (2.3) for 
(7* yield 

{{DUfaDU) = {{DUf) (aDU) = {DUfa^{DU) = a*. (3.11) 
Hence passing to the antisymmetric part it follows that 

/ / -Ai 2(-D;7) -Ai 3(-D(7) \ \ 

(ay = liDUfa^DU) = (c(DUfJDU) = ( c a^^du) o -a^.^du) ) , (3.12) 

\ \Ai,3iDU) A2,3(W) // 

where Aij{DU) := diUid2Uj — diUjd2Ui. Therefore, since cr* is partially isotropic, we obtain 
the following formula for the effective coefficient c*, 

c, = (cAi,2(i^C/)>. (3.13) 
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Using the Cauchy-Schwartz inequality we have 

\\\d2Ui\J \\d1U2\J / (314) 
<c{\d,m\^ + \^2U,\^'y^\\^^U2\^ + \d2U2\^y^\ 
On the other hand (3.11) also imphes for i = 1, 2, 

cr*ei ■ Ci = (aVui ■ Vui) > a{\diUif + \d2Uif). (3.15) 
Combining (3.14) and (3.15) gives the desired bound (3.3). 

Derivation of arbitrarily large coefficients c*; Let 9, k be two positive numbers, let be 
the vectors defined by 

7^' 7^)' 

and let cr^^, cr^, be the (transversely isotropic) phases defined by 

0\ 

al^:=l k9-^ , := /, := 2 1 + J. (3.17) 

V 01; 

Consider the rank- two laminate mixing in the direction C,g the phase ^, with volume fraction 
1 — 6, and the rank-one laminate, with volume fraction 9, composed of the mixture in the 
direction of the phases cr^ and cr^ with volume fraction |. The two-scale conductivity ae^^ is 
defined by 

<^eAy, ^) - XeiCe ■ v) + (l - XeHl ■ v)) {xi^ ■z)a^+{l- x{e ■ z)) a^) , (3.18) 

where y — ^1 z — ^ are the ordered fast variables, xe is the 1-periodic function which agrees 
with the characteristic function of [0, 1 — ^] in [0, 1], and x is the 1-periodic function which 

agrees with the characteristic function of [0, |] in [0, 1]. By [27, 12] the local electric field Eq^^ 
associated with the conductivity (Tth,K has the same laminate structure as (3.18), and thus can 
be written as 

EoAy. ^) Xeiil ■ y) + (1 - xeiil ■ v)) {xie ■ z) + (l - xH' ■ z)) El,) . (3.19) 
The constant matrices E] Eq ^, E^ ^ are the solutions of the linear system 

f (1 - ^) El, + I {El, + El,) = / average-value 

El^ — El^ = i'^ ®ri2 jump of the curl at the scale £^ 

El, - \ {El^ + El^) = (g) jump of the curl at the scale e (3.20) 

{(r'^Eg^ — (T^El^Y^'^ = jump of the div at the scale 

. [(^IkEI, - \ {(y^El^ + a'^El^)fil = jump of the curl at the scale e. 

We refer to [12] for more details. Similarly to (3.11) and taking into account the two-scale 
structure (3.18) the effective conductivity , is given by 

4. = (1 - ^) {El,)^'yl.El, + I [{El,)^a'El, + {El,)^a'El,] . (3.21) 



Taking into account the values (3.17) of the matrix conductivities we deduce that 

(4J^ = ^ {EIJJEI^. (3.22) 

Using Maple to compute explicitly the solutions Eq ^, Eq^ of the linear system (3.20), we 
get the following asymptotics as 6' ^ 0, 

(4«)^ = { ^1,2 (i^U J + 0{e) = -^J + 0(6) (3.23) 
Therefore, cr^^ is asymptotically partially isotropic, and the effective coefficient 

ct^:^^-A,,,{El^)^-^ + 0{9), (3.24) 



is both negative and arbitrarily large when k is arbitrarily large. Moreover, the (2 x 2) deter- 
minant A of the electric field satisfies 

A,,2 (^U = -Ai,2 + 0(1) = + 0(1), (3.25) 

and thus for large k has not the same sign throughout the material, contrary to the columnar 
case. 

On the other hand, if we replace in (3.17) the matrix by 

/2 -1 0\ 

(7^ := 1 2 , (3.26) 

\o y 



then the previous procedure leads us to the asymptotics 

(^;, = ? ^1,2 {ElJ J + 0(6) = ^J + 0{9). (3.27) 



Hence, the effective conductivity al ^ is still asymptotically partially isotropic, and the effective 
coefficient 

cl^=^-A,,,{El^)^-^ + 0{e), (3.28) 
is arbitrarily large but positive. As before, the minor A 1^2 of the electric field satisfies 

Ai,2 (4J - -Ai,2 + 0{l) = ^ + 0(1), (3.29) 

and does not have the same sign throughout the material when k is large. 

4 Hashin-Shtrikman type bounds under geometric isotropy 

4.1 F-tensors, F-operator, and geometric isotropy 

For given a, ^ > 0, consider a periodic two-phase composite with local conductivity 

(^{y) = xi{y) (^i + X2{y) 02 e M«(a, /3; Y), (4.1) 
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where Xi is the characteristic function of the phase i with volume fraction /j, for i = 1,2. 
Denote by a* its effective conductivity. Following [21] (see also [26], Chapter 19) there exists 
an effective tensor associated with the conductivity a, defined by (e e Lj(y)^ is the electric 
field and j e Lj(F)^ is the current field) 

P(i) = -nP(e), where j^ae and P{g) -.^ {xi{9 - (g))). (4.2) 

In some sense P is the projection on phase 1 of the fluctuating component of the field. Also 
we have for the adjoint problem 

P{j')^-Y,^P{e'), where / = aV. (4.3) 

Recall the relation (2.8) and the definition (2.7) of the tensor L{y) which enters it. A similar 
computation based on the formulas (4.2) and (4.3) lead us to an effective tensor associated 
with the tensor L{y) of (2.7), and defined by 

where similarly to (2.10), 

Now, we will derive a Hashin-Shtrikman type variational inequality associated with the 
variational principle (2.9). To this end, let us consider for a given reference tensor Lq, the 
nonlocal operator F defined for periodic vector- valued functions A,B& Lj{Y)^, by 

B^TA if TiB ^ B and ri{A- LoB) ^ A- LoB, (4.6) 

where Fi represents the projection on the space of fields which satisfy the same differential 
constraints as E E L|(F)® in (2.8). Since E is composed by a divergence free field js G (F)^, 
and a curl free field £ Lj{Y)^, the operator Fi in Fourier space is given by 



Fi(A;)=Fi(0= (^^ where ^ := -, for A; e {0}. (4.7) 

Under the conditions Li > Lq > 0, for i = 1,2, the Hashin-Shtrikman type variational inequality 
associated with the variational principle (2.9) is given by the formula (13.30) of [26], which reads 
as 

(L,-Lo)~'(F):(F)<([F + (L,-Lo)-']F:F>, for any F e Lj(y)^ (4.8) 
Following the computations of [26] (Section 23.6) this inequality implies the bound 

-\ -1 



tV E \xiik)\"m 

■'^•'^ fcez3\{o} 



(4.9) 



which also holds for the enlarged inequalities Li > Lq > 0, for i = 1,2. Note that by virtue 
of the Plancherel equality the Fourier coefficients Xi{k) of the characteristic function xi satisfy 
the equality 

tV E \uk)\" = j^{{xi-fir) = i. (4.10) 

■'^■'^ fe€Z3\{0} ■'^■^^ 
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So, the series in (4.9) can be regarded of an average of the operator F. 

Finally, consider the case of a two-phase random composite. According to [34] (see also 
[26], Section 15.6) the composite is said to have a geometric isotropy if all correlation functions 
associated with the geometry represented by the characteristic function Xi ^ire invariant by 
rotation (or reflection). Then, under geometric isotropy the average series of (4.9) reduces to 
an average of F over all directions of the unit sphere S'^. Therefore, we get the bound (see [26], 
Section 23.6) 

y* + i^o>7^. (4.11) 
4.2 Hashin-Shtrikman type bounds 

Consider a periodic two-phase composite with non-symmetric positive definite conductivities 



(7i := Q tti , with bi > 62, (4.12) 




of respective volume fractions /j, for i = 1, 2. The effective conductivity cr* of the composite is 
assumed to be transversely isotropic, i.e. 



(4.13) 




Let g : (0, 00) — >■ IR be the function defined by 



9ir)---7;l ...r, ■-712. ^^ g (0.1). forr>0. (4.14) 



1 r cos^^sing 

2 Jo cos2 e + r-1 sin^ 9 

Consider the coefficients a±, tf, sf , ay, by, cy defined by 



aiC2 - a2Ci ± J aia2 ((ai - 02)^ (ci - C2Y) 

a± := ^- , (4.15) 

Oi - 02 



t( ^ ■ \ f( , • ^ , /1/2 (ai + «ci - a2 - ic2)^ , . 

ay + icy := -/2(ai + ^ ci) - /i(a2 + 102) + j— . ^ , , . — ^ (4.17) 

(/i(ai + ici) + /2(a2 + IC2) - a^-ic^) 



/1/2 {hi - b 



\2 



JlOl + /2O2 - 0* 
Then, we have the following result: 

Theorem 4.1. Assume that the composite is geometrically isotropic. Then, in view of def- 
initions (4-15)-(4-18) the coefficients a*, c* of the effective conductivity cr* (4-13) satisfy the 
Hashin-Shtrikman type bounds 

a^y + (cy + a±y-^<0, (4.19) 
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c 




while the points (ai,Ci), (02,02) solve the equations 



a 



a^ + (c-a±)^- — = 0. (4.20) 
^1 

Moreover, the coefficient 6* satisfies the hounds 

111, . ^ 

^ + ^ > -, — > 0. 4.21 

Remark 4.2. In the ay-cy plane the bounds (4.19) correspond to the intersection of two disks 
parametrized by a±, which are tangent to the cy-axis. Due to definition (4.17) these bounds 
remain the same if we replace Ci, C2, — cy by ci + Cq, C2 + Cq, — cy + Cq. This refiects the fact 
if we add a antisymmetric matrix to the local conductivity cr, then the same antisymmetric 
matrix is added to cr*. Also note that if ci = C2 = c* = c, then cy = — c. 

Remark 4.3. With the change cy to — cy, the circle 0±^ satisfying the equality in (4.19) is 
the same as the circle 0± of equation (4.20) passing through the points (oj, Ci), i = 1, 2, when 
sf = tf. Moreover, since g{r) G (0, 1) for r > 0, we have < sf < tf in (4.16). This implies 
that the radius (2s^)~^ of 0±^ is larger than the radius {2tf)~'^ of 0±- The circles Qi^ 
0± are also tangent at the same point of the c-axis. The geometrical picture is given by 
Figure 1 in the a-c plane. 
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Remark 4.4. The inequalities (4.19), (4.21) do not allow us to show that c* tends to zero 
when ci and C2 approach zero, while keeping oi, 02, &i, and &2 fixed. To this end, we only have 
the bound (3.9) which reads as 



1 2 — Ci — C2 1 ^ 



min (ai, 02) 



Cl - C2 . 



(4.22) 



Proof of Theorem 4.1. The proof is divided in four steps. In the first step we determine 
a suitable reference tensor Lq. In the second step we compute the tensor r(^) involved in the 
y-tensor approach. In the third step we compute the average (r)52. The fourth step is devoted 
to the derivation of the bounds. 

First step : Determination of Lq. 

Similarly to (2.7), let Lj, for i = 1, 2, *, be the symmetric tensor defined by 



.<yt (^f 



^f-^f(^f)-V 



/i 











1 

















1 

bi 





_ Ci 
Ui 





Q 








ai 






\o 













- 



ai + 






ai 







o\ 







ai + — 

hj 



Now, let Lq be the symmetric tensor defined by 



^0 •= 



Cl C2 



/ti 



t4 

-ts 

t2 

\0 



t^ 0\ 

-^2 



t3 

t3 

ts/ 



where Cj e 



-Co 



The condition Lq > is equivalent to 

tx > 0, ^4 > 0, ^3 > 0, ^5 > 0, txt^ > tl 
We also need Li> Lq, which is equivalent to 

/I 



> tl, det 



ai 

Ci 



-tl 



Ci \ 

12 

ai 



\ ^2 C^i H ^3 

\ai ai 



> 0, for i = 1, 2, 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



— > ^4, bi > ^5, for i = 1,2. 

Oi 



Prom now on, assume that bi > 62, and set 

t4 := 



1 



62, 



(4.27) 



(4.28) 
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in order to make the inequalities (4.27) as sharp as possible. 

On the other hand, the inequalities (4.26) show that the points (ai,Ci) and (02,02) belong 
to the disk in the a-c plane, 

ti (a^ + c^) + (4 - tih - 1) a - 2*2 c + < 0, (4.29) 

which lies in the half-plane a > 0. To make these bounds as tight as possible, we consider 
the two circles which are tangent to the c-axis, and which pass through the points (oi, Ci) and 
(02,02). This requires 

h = ^, (4.30) 

and the two circle equations 

tjiaf + Ci) -tia,,-2tit2Ci + tl = 0, fori = 1,2, (4.31) 
which can be written as 

t? 02(0? + c?) - ti 01O2 - 2tit2 02C1 + i| 02 = 

(4.32) 

tf Oi(o| + c|) — ti O1O2 — 2tit2 O1C2 + t| Oi = 0. 

Subtracting and dividing by tf, we get that a := ^2/^1 solves 

(oi - 02) - 2 (01C2 - 02C1) a + oi(o2 + cl) — 02(0^ + cf) = 0, (4.33) 

the discriminant of which is 

(01C2 -02Ci)^ + (oi - 02) (02(01 + Ci) -01(02 + 02)) = 01O2 ((oi -02)^ + (ci -02)^) > 0. (4.34) 

Hence, equation (4.33) has two real solutions (one for each circle) a± which are given by (4.15). 
Moreover, putting ^2 = 0;^! in (4-31) we obtain that 

^^^^-r^ for z = 1,2, (4.35) 

which implies that the (2 x 2) matrix in (4.26) is non-negative. Therefore, the choice of the 
coefficients ^1,^2, ^3 given by 

^1 = ^ 7-^ ^2 = 0;^!, ^3 = 0:^^1, for q; = q;±, (4.36) 

of + (ci - ay 

combined with (4.28), implies the desired inequalities Li > Lq, for i = 1,2. Making this choice 
in (4.29) the points (ai,Ci) and (02,02) belong to the two circles of equation (4.20) which are 
tangent to the line o = 0. 

Second step : Computation of r(^). 

Let ^ e 5"^, ^ = (sin 9 cos </?, sin 9 sin (/?, cos 9). By virtue of Section 4.1 the tensor r(^) is defined 
from the tensor Lq (4.24), by 

(t) ^ ^^^^ ) ' ^ ^^-^^^ 

if and only if 
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(4.38) 



and 



e<H)e 



A2 



Ci C2 



Bi 
B2 



(4.39) 



By (4.38) we have B^^ — 0, and = ^ (g) for some vector r]. Prom (4.39) it follows that 



(4.40) 



(4.41) 



Ali - BlC2i - BlCli = Ali - BlC2i - (Cse ■ = 0, 
- C,B^ - C2 (e ® 77) ^i® {Ali - Bleu - B^CU) 

= e ® {AJ^ - BfCU - {C2^ ■ V) 
= ^ ® (Af e - BfCf^) since C^^ = - C2, 
= e ® where A; := Aj^ - BjCj^. 
Noting that Cf^C2 is antisymmetric, this implies that 

= = (<^r'A)^e-(^®e) (c^r'c^2)'^e-(^®e) (c'r')^e = (Q-'A)'^e-(Q-'e-e) ^, (4.42) 

so we have 



{C^'A,)U 



(4.43) 



Moreover, replacing Bi given by (4.41) in (4.40) and using that (C^ ^)^C2 is antisymmetric, we 
get that 



O^AU- {C^'A^fC2C + iC{'C2fC2C + (A; ® {C{TC2C - iv ^ cjc 



AU - {C^'A^rC2C - [{Cs - CjC^'C2) ^■^]v, 



hence 



V 



A^ - {c^cr'A.n 



where D := C3 - CjC^^C2. 



Again using (4.41) combined with (4.43) and (4.45) we deduce that 



(4.44) 



(4.45) 



Bo 



- 1. \ C2 Ml. 



(4.46) 



Hence, from definition (4.37) it follows that 



/^-i ^ c^'c2 (e ® ^^2"^^^-^ (e ® ^^r' ^^r'c^2'^(e ® o \ 

i 



r(0 = 



\ 



which is a symmetric matrix since Cf = Ci and = — C2. 
Third step : Computation of ((r)52)~^. 

Note that the computation of r(^) can be carried out if the matrix D of (4.45) 



(4.47) 



D 



di 











d. 











d2 



h- 






ti 











\ 



^3-^ 





(4.48) 
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is positive definite, i.e. tit^ > t\ and ^5 > 0. Let us assume these conditions for the moment. 
We shall be able pass to the hmit as di — > in the expression of ((r)52)~^. Set 



Q 



'qi 0' 
gi 















Pi 












i) 




R 




52 





' 

n I := - C^^QC^K 

r2, 



(4.49) 



(4.50) 



By definition (1.6) we have 



Pi 



1 r sin^e 



d2 cos^ 9 



dd — 00, 



P2 



^2 

S3 



dl + {d2 — di) ^3 / g2 di-^o d2 t5 



(4.51) 
(4.52) 



Moreover, the matrix 



ti 









qi 





ti 











^4/ 



(^ cr 



(4.53) 



52 



has the property that its trace is 1. This combined with definitions (1.6) and (4.14) yields 



qi 



(4.54) 



which also implies that 



1 _ ^ = ± 

1 92 1 



(l + ^(il/t4)) 



(4.55) 



14 ^4 t4 



On the other hand, by definition (4.49) we have 

a2 



^ 2 
|pi 

\ oy 



- 



ti 
\ 



-Pi 0^ 

ti 



oy 



(4.56) 
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Then, putting (4.50) and (4.56) in the S'^-average of (4.47) we get that 







t2 







tl 



t2 












u 


t2 




r2 








Pi 















t2 A 



which gives 



/I 

r-i 






t2 

tin 





1 

n 



t2 





1 

r2 



t2 

tin 





f2 



tin 




-+,2 

Pi tf n 








pi 



h_ 
tin 







''2 















1 

Pi ri 










P2J 










1 

P2/ 



(4.57) 



(4.58) 



Therefore, passing to the hmit as di — > 0, or equivalently ^3 — )■ tl/ti, (4.51) and (4.52) imply 
that 



/I 











^2 


°1 


ri 








tin 






1 




t2 
























ri 


1 


tin 



























r2 











^2 





''2 










tin 










t2 
















tin 








tjn 




\ 

















(4.59) 



Fourth step : Derivation of the bounds. 

On the one hand, the Appendix of [21] (see also formula (19.3) of [26]) yields the following 
formula for the F-tensor defined by (4.2) 



= -/2 (Tl - /l (T2 + /1/2 ((Tl - (T2) (/i (Tl + /2 (T2 - (T*) ((Ti - (T2) 

Note that, due to the transverse isotropy of (jj, for i = 1, 2, *, we have 

^ay —cy 

Y^= \ Cy ay 




(4.60) 



(4.61) 
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This relation also separates into blocks, so that we obtain for the 33 entry of the relation 
(4.18). Moreover, making the correspondence 



a —c 



i — > a + ic, 



we deduce from the first (2 x 2) block of (4.60) the relation (4.17). 
On the other hand, by (4.61) the formula (4.5) for V* reads as 



(4.62) 



1 











Cy 


o\ 


ay 








ay 







1 





Cy 










ay 


1 


ay 






































Cy 





Cy 

ay + — 










ay 




ay 






Cy 











Cy 

ay + — 





ay 








ay 



















by) 



(4.63) 



Then, the bound (4.11) applied with the formulas (4.63) for (4.24) for Lq and (4.59) for 
((r)52)"\ combined with (4.28), (4.54), and (4.55), implies that 



1 1 1 

+ ^ > 



by bi r2 &i(l - 5'(Mi)) ' 



by > 0, 



(4.64) 



and 



det 



ay 



Cy \ 

\-S2 ^ 

ay 



Cy Cy 
h S2 ay -\ S3 

,ay ay 



>0, 



(4.65) 



where 



Sl - 



S2 



S3 



1 

n 

t2 



tin 



-to 



l + (7(Mi) 

2 

l + ^(Mi) 



- 1 ) > (since < g{biti) < 1) 



(4.66) 



tin 



-t. 



Due to (4.30) we have S1S3 = s^. Therefore, similarly to (4.26) and (4.29) the inequality (4.65) 
can be written as 



ay + [ Cy + 



S2 
Sl 



ay 

Sl 



< 0. 



(4.67) 



Finally, taking into account (4.16), (4.36), (4.66) the inequalities (4.67) and (4.64) correspond 
respectively to the desired bounds (4.19), (4.21). Theorem 4.1 is proved. □ 
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